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Large Gravitational Wave Background Signals in Electroweak Baryogenesis Scenarios 
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The bubble wall velocity in an electroweak first order phase transition is a key quantity both for 
electroweak baryogenesis and for the production of a stochastic background of gravitational waves 
that may be probed in the future through gravitational wave experiments like LISA or BBO. We 
show that, contrary to the conclusion drawn from previous studies, it is actually possible to generate 
a potentially large gravitational wave signal while satisfying the requirements for viable electroweak 
baryogenesis, once the effects of the hydrodynamics of bubble growth are taken into account. Then, 
the observation of a large gravitational wave background from the electroweak phase transition 
would not necessarily rule out electroweak baryogenesis as the mechanism having generated the 
observed baryon asymmetry of the universe. 

PACS numbers: ll.30.Fs, 47.35.Bb 



I. INTRODUCTION 

The study of bubble growth in cosmological first or- 
der phase transitions is very relevant for many phenom- 
ena having possibly occured in the early universe such 
as electroweak baryogenesis [l|, [1| or the production of 
a stochastic background of gravitational waves 
first order phase transition proceeds by bubble nucleation 
and expansion, and the velocity of the expanding bubble 
walls plays an essential role both in electroweak baryo- 
genesis and in gravitational wave production, since the 
efficiency of both processes strongly depends on its value. 

Treatments of the bubble wall velocity generally as- 
sume that friction from the plasma balances the initial 
pressure difference that drives the bubble expansion, so 
that the wall reaches a constant speed after a short pe- 
riod of acceleration. Assuming that the free energy of 
the Higgs field is released into the plasma, a hydrody- 
namic treatment of the plasma can be used to determine 
the fluid motion 0| , but this approach leaves the wall 
velocity as a free parameter as long as the microscopic 
mechanism of friction is unknown. Ultimately, the wall 
velocity can be fixed using the equation of motion for 
the Higgs field, that takes into account the friction of the 
plasma 8., .9]. 

However, successful electroweak baryogenesis and siz- 
able gravitational wave production require very different 
wall velocities. The electroweak baryogenesis mechanism 
is based on the interaction between the expanding bubble 
wall and the plasma in front of it, leading to a CP asym- 
metric reflection on the wall of certain particle species, 
and the subsequent diffusion of these particle asymme- 
tries into the plasma in front of the bubble wall [2|] , where 
sphalerons are active and capable of converting the CP 
asymmetry into a net baryon number. Then the gener- 
ated baryon number is carried into the broken phase as 
the wall passes by (where it stays frozen if the sphaleron 
processes are sufficiently suppressed in the broken phase). 
In order for this whole mechanism to be effective the dif- 
fusion timescale has to be smaller than the time the wall 
takes to sweep through the plasma just in front (other- 



wise diffusion is ineffective and the generation of baryon 
number is strongly supressed), and this puts an upper 
bound on the relative velocity between the wall and the 
plasma in front V < D/L^ ~ 0.15 — 0.3 (being D a cer- 
tain diffusion constant and L„ the wall thickness) 
Moreover, it is generally stated that subsonic wall veloc- 
ities are always needed, because effective diffusion can- 
not take place for supersonic walls. On the other hand, 
fast moving walls are essential for the production of a 
sizable amount of gravitational radiation in bubble colli- 
sions d, [Tol - [l^ , or turbulence in the plasma [H, [l3] . 
In particular, for bubble collisions the gravitational wave 
amplitude is roughly proportional to [l^, Ull , with an 
extra implicit T4;-dependence through the efficiency co- 
efficient K for transforming the available energy from the 
phase transition into plasma bulk motion (which is in 
turn responsible for the generation of the gravitational 
wave background during the bubble collisions 0, [lH), 
and this dependence further suppresses the gravitational 
wave signal for small wall velocities (this suppression also 
affects the gravitational wave signal generated from tur- 
bulence). All this has established the common lore that 
both phenomena cannot happen in the same scenario. 

In the Standard Model, and for values of the Higgs 
mass above the LEP bound Mh > 114.4 GeV [2fl], the 
electroweak phase transition is found not to be of first or- 
der, but rather a smooth cross-over [2l|. However, there 
are many possible theories beyond the Standard Model in 
which the electroweak phase transition may naturally be 
of first order, such as extensions of the MSSM [H-Q (in 
the MSSM itself, the region of parameter space where a 
first order phase transition leading to electroweak baryo- 
genesis is achieved is currently very tightly constrained 
|25|). Two-Higgs-Doublet model s [26l-[29| . singlet field ex- 
tensions of the Standard Mod el I30l434| . composite Higgs 
models and others. Also, in |35| the electroweak phase 
transition was studied for the Standard Model considered 
as an effective theory with a low cut-off, finding that the 
inclusion of higher dimensional operators in the Higgs po- 
tential may give rise to a rather strong first order phase 
transition. As it has been discussed above, for these or 
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any other model leading to a first order electroweak phase 
transition, the value of the velocity of the expanding bub- 
bles is a key parameter for the study of both electroweak 
baryogenesis and gravitational wave production at the 
phase transition. For the case of the MSSM the value 
of Vw was found to be quite small over all the available 
parameter space [36|, but it is expected that it may be 



much larger in many of the models discussed above, since 
the electroweak phase transition in those cases is much 
more strongly first order, since the wall velocity increases 
with the strength of the phase transition [g] . 

Here we will show that in contrast to the case of grav- 
itational wave production, where the relevant velocity is 
indeed the speed of the wall V^, in electroweak baryo- 
genesis the relevant velocity (being the relative velocity 
between the bubble wall and the plasma just in front 
the wall v+) is in general lower than Vw, and this ef- 
fect becomes more important as the phase transition gets 
stronger. Then, it is possible to have a sizable gravita- 
tional wave production (through a relatively large Vw, a 
natural possibility in many of the scenarios beyond the 
Standard Model mentioned above, where the electroweak 
phase transition can be strongly first order) while the 
electroweak baryogenesis mechanism is still effective (w+ 
is sufficiently low). 

This article will be organized as follows: In section 2 
a summary of the hydrodynamic analysis of bubble ex- 
pansion is given. Then, in section 3 the use of v+ instead 
of Vw as relevant velocity for electroweak baryogenesis is 
motivated, and in section 4 it is shown that v+ < Vw (and 
v+ <C Vw for strong transitions). In section 5 the gravita- 
tional wave amplitude in viable electroweak baryogenesis 
scenarios is obtained, and we conclude in section 6. 



II. HYDRODYNAMIC RELATIONS. 

The hydrodynamic analysis of the system consistin g on 
the Bubble wall and the plasma surrounding it ^, [tI Il5j| 
rests on two basic assumptions: energy-momentum con- 
servation in the system and local thermal equilibrium 
(LTE) in the plasma. This last assumption is reasonable 
since local equilibration due to the interactions between 
the particle species in the thermal plasma is much faster 
than the actual macroscopic movement of the plasma. 
Local thermal equilibrium is also crucial for the consis- 
tency of the fluid approximation applied in electroweak 
baryogenesis computations (see 0). 

The energy-momentum tensor of the Higgs field 4> is 
given by 



-dpC^dPct^^Voic^) 



(1) 



where Vo((/>) is the renormalized vacuum potential. If 
the plasma is locally in equilibrium its energy-momentum 
tensor can be parametrized as 



rpplasma 

fJ.L/ 



(2) 



where w and p are the plasma enthalpy and pressure, 
respectively. The quantity is the four-velocity field 
of the plasma, related to the three-velocity v by = 
(7,7v). The enthalpy w, the entropy density a and the 
energy density e are defined by 
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where T is the temperature of the plasma. Conservation 
of energy- momentum is given by d^Tj^^ + jriasma _ 
Since we are interested in a system where the bubble 
expands at a constant speed, energy-momentum conser- 
vation reads in the wall frame (assuming no time depen- 
dence, and with the wall and fluid velocities aligned in the 
z direction) dzT''-^ = dzT^° = 0. Integrating these equa- 
tions across the wall and denoting the phases in front and 
behind by subscripts -|- (symmetric phase) and — (bro- 
ken phase) one obtains the matching equations in the 
wall frame: 



w+v'^j'^ + p+ ~ w^v'^^'^ + P- 



(4) 



From these equations we can obtain the relations Q 



■P+ 



P- 



(5) 



In order to proceed further, one needs to know the equa- 
tion of state (EoS) for the plasma. A parametrization 
that accounts for deviations from the so-called bag EoS 
(usually used in hydrodynamical studies of phase transi- 
tions in the early universe d, 0|) is (see [Tsl ) 



(6) 



P- = \^-^- + ^ = «+^+ + 
e_ = a-Tt — e = a+T^ i 



where we have defined a+ = 



and r 
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. The 



quantity a+ is approximately the ratio of vacuum energy 
difference to thermal energy in front of the wall, and char- 
acterizes the strength of the phase transition (the larger 
Q!+ the stronger the transition) , and a± are related to the 
number of relativistic d.o.f in the symmetric and broken 
phases. Using ([6]) and ^ we can write the relations ([5]) 
as 



1 - (1 - 3a+)r 
3-3(H-a+)r 



3 -I- (1 — 3a+)r 
V- ~ 1 + 3(1 + a+)r 



(8) 
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The two equations 
1 

v+ = 



l + a+ 



can be combined to give 
1 



2 ^ 



(9) 



so that there are two branches of solutions, correspond- 
ing to the ± signs in ([9]). 

In a concrete model where p = —T (the free energy 
or finite-temperature effective potential) is known, the 
thermodynamic potentials can be calculated in the two 
phases and the temperature T/v at which the phase tran- 
sition happens is determined using standard techniques 
[l6]. Still, there are three unknown quantities (T-, w+ 
and u_ ) and two equations ([5]) , so up to this point all hy- 
drodynamically viable solutions are parametrized by one 
parameter, usually chosen to be the wall velocity Yw 

Next, we briefly review the solutions of the plasma 
velocity v \v\ . Applying energy- momentum conser- 
vation in the plasma (}Mypiasma _ q (^^^ from the wall 

j^^t gives a constant background that plays no role 
in the energy-momentum conservation), we obtain a set 
of hydrodynamic equations. Since there is no intrinsic 
macroscopic length scale in the system, solutions to these 
equations are self-similar and only depend on ^ = r/i, 
with r being the radial coordinate of the bubble and t 
the time since nucleation. The plasma then fulfills the 
equations 



2- = 7^(1 - vi) 
w 



^-1 
ci 



(10) 



where Cs = 1/a/3 denotes the velocity of sound in the 
plasma and iJ,{^,v) is the Lorentz-transformed fluid ve- 
locity ^(^, v) = jE^- Generally, there are three different 
types of solutions to pH]) with the boundary conditions 
© [3 [13 ■ detonations, deflagrations and hybrid so- 
lutions. In detonations the bubble wall expands at su- 
personic velocities and the vacuum energy of the Higgs 
leads to a rarefaction wave behind the wall, while the 
plasma in front is at rest. In this case, the wall ve- 
locity is Vjj] — v-f- > v^, and therefore detonations are 
identified with the -I- branch of solutions in ([S]). In de- 
flagrations, the plasma is mostly affected by reflection 
of particles at the bubble wall and a compression wave 
builds up in front of the wall while the plasma behind 
is at rest. In this case, the wall velocity is identified 
with Vw = V- > v^, corresponding to the — branch 
of solutions in ([9]). While "pure" deflagrations are sub- 
sonic, the hybrid case occurs for supersonic deflagrations 
where both effects (compression and rarefaction wave) 
are present. From now on we will focus on deflagrations 
and hybrids since for them > «+. Also, in this case 



T+ > Tn, and so a+ < aN 



(being a ^ the actual 



III. v+ VS Vn, IN ELECTROWEAK 
BARYOGENESIS. 



Consider the evolution of a plasma volume element as 
the compression front and bubble wall reach it and pass 
by, in the case of a deflagration (Figure [Ij. The ele- 
ment is initially at rest (its position r — cte), and as it 
enters the compression wave it acquires a velocity, that 
grows until reaching v{Vw) = close to the bub- 

ble wall. Then, due to the compression wave, the rela- 
tive velocity between the wall and the volume element 
just in front of it is w+ independently from the details 
of the electroweak baryogenesis mechanism close to the 
wall, just relying on LTE applying all over the compres- 
sion wave. In the electroweak baryogenesis analysis the 
relevant velocity is the perturbation V with respect to 
the background velocity of the plasma w+. However, in 
[1] (and subsequent electroweak baryogenesis studies) the 
background velocity of the plasma was taken to be V^, 
neglecting hydrodynamics. We find that in the presence 
of a compression wave and in the wall reference frame, 
the velocity of the plasma should he V = — ?;+ -I- V in- 
stead oi V — —Vu, + V, and so the dependence of the 
electroweak baryogenesis mechanism on the wall velocity 
Vu, extracted from is actually a dependence on v+. 
The same argument applies to the upper bound on Vy^ 
from [2I from the requirement of effective diffusion of the 
CP asymmetric current in front of the wall. The bound 
should now be regarded as w+ < D/L^ ^ 0.15 — 0.3 
(the wall thickness ^ (15 — 40)/r and the diffusion 
constant Z? - 5/r [2]). 




Electroweak Baryogenesis 
Happens Close to the Wall 



measure of the strenght of the phase transition). 



FIG. 1: Movement of plasma volume elements with time. 



IV. DEPARTURE OF v+ FROM K,. 

In the presence of a compression wave (deflagrations 
and hybrid solutions) it is seen from dH) and (fTO|) that 
Vw > vj^. For weak first order phase transitions (aAr <C 
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1) one has ~ 14; from ([9]) and the effect of considering 
Vj^ instead of as the relevant velocity for electroweak 
baryogenesis is small. However, as a a? gets larger (al- 
ways keeping a+ < 1/3 Ts']) and the phase transition 
gets stronger, w+ progressively departs from Vw, even- 
tually reaching v+ ^ Vw for very strong phase transi- 
tions. This is shown in Figure [21 where v+{Vw) is plot- 
ted for increasing values of aiq- Then, for rather strong 
phase transitions, v+ can be kept small enough to satisfy 
the diffusion upper bound for electroweak baryogenesis 
v+ < 0.15 — 0.3 with a rather large 14, (T4 ~ for 
deflagrations or even > Cs for hybrid solutions). 



V, 0.2 - 




FIG. 2: Relation v+{aM ,Vw 



In Figure [3] the region in the parameter space of ajv 
and Vw compatible with electroweak baryogenesis (for 
various values of the upper bound on vj^) is shown. For 
large aM and small 14>, a+ > i and no solutions exist 
[l5]. Also, for small aN there is a maximum value Vw 
can take with positive plasma friction [l8l |. and ajsi has 
to be larger than a critical value ac for bubble expansion 
to be possible (etc — 0.05 for a_/a+ = 0.85) lisj. 






= 0.15 





= 0.2 





= 0.3 



No Solutions 



«N 0.5 




FIG. 3: Regions in the (T4j,ajv) plane compatible with v+ < 
0.15, i;+ < 0.2 and i;+ < 0.3 (region above each line). 



V. HOW BIG CAN THE GRAVITATIONAL 
WAVE SIGNAL BE? 

Here we will concentrate on production of gravitational 
waves through bubble collisions during an electroweak 
first order phase transition (the analysis can be extended 
to the case of turbulence, with similar conclusions). The 
amplitude and peak frequency of the generated stochastic 
spectrum are [ll| (see also [sl. [To|): 



i^Gwh (fcoii) ^ 10" 



100 

9* 



ajv 



l + Qjv/ 0A2 + V3 



fcoli 



10 



Vloo; 



T 



100/ VlOOGeVy H,1.8 + V;i 



1.2 



(11) 



(12) 



and the spectrum grows as for frequencies smaller 
than J^oi; [HI, [l^l and falls off as for large frequen- 
cies [ul. Typically P/H^ - 100 [H and for the elec- 
troweak phase transition T ~ 100 GeV and ~ 100. 
Then, if viable electroweak baryogenesis would require 
Vw < 0.15 — 0.3, this would suppress the gravitational 
wave amplitude both explicitly and through the effi- 
ciency coefficient K{aN,Vw) dependence even for rather 
strong phase transitions {un ^ 0(0.1 — 0.2)), making the 
gravitational wave signal in scenarios where electroweak 
baryogenesis is possible undetectable both at LISA and 
BBO (blue lines in Figure S]). Moreover, for very strong 
transitions {aN C(l)) it would be impossible to satisfy 
the bound < 0.15 — 0.3 and electroweak baryogenesis 
would simply not be possible. 



le-lO 




FIG. 4: flaw h'^ if) for various values of qjv and Vw = 0.15 
(blue lines) or i;+ — 0.15 (red lines). 



However, since the electroweak baryogenesis bound ac- 
tually applies to w+, the gravitational wave signal am- 
plitude for a scenario where the bound is satisfied gets 
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enhanced by 2-3 orders of magnitude with respect to the 
previous situation (for a transition of the same strength) , 
and it is then possible to achieve electroweak baryogene- 
sis in scenarios where the gravitational wave signal may 
be detectable at BBO for moderately strong phase tran- 
sitions (solid red lines in Figure |4]) or even marginally at 
LISA for extremely strong transitions (dashed red lines 
in Figure 13]). 

VI. CONCLUSIONS. 

If the electroweak phase transition is of first order, it 
will proceed by bubble nucleation and expansion. In the 
presence of a compression wave in front of the expanding 
bubble walls (which always occurs if the bubbles expand 
subsonically and can also occur under certain conditions 
if the bubbles expand supersonically), the relative ve- 
locity between the wall and the plasma in front («+) 
is smaller than Vw While this effect is small for weak 
phase transitions (for which one gets v+ ~ Ku), it be- 
comes important as the phase transition gets stronger, 
and for rather strong transitions one has <C • Since 
the background velocity of the plasma relevant for the 
electroweak baryogenesis process is whereas the rele- 
vant velocity for gravitational wave production through 
bubble collisions is Vw (see however this opens the 

possibility of achieving viable electroweak baryogenesis 
(satisfying the diffusion bound v+ < 0.15 — 0.3) and a 



sizable gravitational wave signal in the same scenario, 
which was previously regarded as not possible due to 
the very different ranges of wall velocities T4, that were 
thought to be required for each the two processes to be 
efficient. We actually find that in electroweak baryogen- 
esis scenarios, the gravitational wave signal can be eas- 
ily detected by BBO for moderately strong phase tran- 
sitions (ttAT ~ 0(0.2 — 0.3)), and is very close to the 
LISA sensitivity curve for extremely strong transitions 
(«jv^O(l)). 

However, whereas a moderately strong first order phase 
transition is a natural possibility in many theories beyond 
the Standard Model (and so a positive signal at BBO 
is plausible in these scenarios), the occurence of such a 
strong transition as to observe a signal with LISA is quite 
unlikely since it corresponds to a very fine-tuned scenario, 
and also in this last case the bubble expansion would have 
to proceed through hybrids (Vw > Cs) and it has been 
argued that hybrids are not in general stable solutions, 
but may evolve into detonations [131 ■ should have all 
these issues in mind when dealing with a concrete model. 
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